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Properties of the Jacobi i^a-function and its derivatives under discrete Fourier transforms are 

investigated, and several interesting results are obtained. The role of modulo A'^ equivalence classes 

J5^ ' in the theory of i9-functions is stressed. An important conjecture is studied. 
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>>■ I. INTRODUCTION 



m 



Methods in mathematical physics usually provide an interface between quite different areas of physics, and it is 
not unusual that such areas advance in parallel, mostly ignoring the each other's steps. This is the case with finite 

m ■ 

^ ' dimensional inner product spaces (hereafter mentioned as the "discrete"), with its leading role in Quantum Mechanics 

00 ; 

^-H (hence Quantum Information Theory) and in finite signal analysis. References [1-3] provide some links between those 

O 

■TJ" theories). 

o ; 

\^ Both quantum mechanics of finite dimensional Hilbert spaces and finite signal analysis rely heavily on the Discrete 

o 

'-^' Mechanics, after the seminal work of Weyl on finite dimensional systems [4], it was Schwinger who observed and 
-I— > 



Fourier Transform (DFT, sometimes mentioned Finite or Fractional Fourier Transform), and, regarding Quantum 



explored the fact that two physical observables whose families of eigenstates are connected via DFT share a maximum 
degree of incompatibility [5] . 

Although, at first glance, a finite system might look much simpler than anything defined on a non enumerable 

^^ ' infinite dimensional Hilbert space (hereafter referred to as the "continuum" ) , there is much more knowledge about 

5t 1 the later than the former. In one phrase, in the continuum we have one, and only one, harmonic oscillator, while 

in the discrete there are a lot of candidates for that role, each one surely with its virtues, but surely no undisputed 

champion. 

The eigenstates associated to the harmonic oscillator, the Gaussian function and the Hermite polynomials, have a 
very distinguishable behavior under the action of the (usual) Fourier transform, so widely known that any comment 
on this regard is completely superfluous. Over such properties rests a huge amount of physical knowledge. On the 
other hand, however, although the Discrete Fourier Transform (DFT) is a well known tool, there is nothing on this 
context which could claim for itself a role analogous to that of the Gaussian function/Fourier transform "duo" . 

A decisive step in an attempt to "regain" , in the discrete, all interpretative power derived from the qualitative 
behavior of the harmonic oscillator eigenfunctions, lost when one leaves the continuum realm, was given in [6], where 
the eigenstates of the DFT are obtained. The purpose of this paper is to further explore this path, showing new results 



which closely parallel those of the continuum. Those results are obtained in a strikingly simple fashion, exploring the 
technique of breaking infinite sums in modulo A^ equivalence classes. Pertinent research on the eigentates of the DFT 
can also be found in [7]. 

A remark must be made about the orthogonality of the DFT's eigenstates. Mehta has conjectured that those states 
are indeed orthogonal, what seems to be most reasonable. One may be led to believe that, just as in the continuum, 
the eigenstates of the DFT may be also (non-degenerate) eigenstates of some other (unitary or self-adjoint) operator, 
and thus orthogonal. Further evidence supporting such conjecture is that the continuous limit of the DFT eigenstates 
recovers, as expected, the Gaussian times the Hermite polynomials. However, as it will be shown, quite surprisingly, 
the conjecture does not hold, giving another fine example of the peculiarities of the finite dimensional context. 

The eigenstates of the DFT are seen to be the Jacobi ??3-function and its derivatives [8]. Interest in Jacobi -d- 
functions, by their own turn, may come from a variety of directions. First, its mathematical interest goes without 
saying (see, for example, [9] and references therein). To cite relatively recent examples in physics, in quantum physics 
it is deeply related to coherent states associated to both circle [10] and finite lattice topology [11]. Its modular 
properties have proven to be of fundamental importance in superstring theory, as it is shown by standard literature 
in this field [12]. 

The basic notation adopted in this article and some preliminary results are presented in the next section. Following, 
orthogonality of the DFT's eigenstates is discussed. Section IV contains the main results, for which a two variable 
generalization is verified in the subsequent section. Further relations among i?3 functions are obtained in section VI, 
which precedes the concluding section. 

II. PRELIMINARY RESULTS 



In reference [6] it is shown that there is a set of functions with the following remarkable property 

.„ N-l 
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fnij) = -7= XI /n(^)exp 
^^ fc=0 

where A'^ is a natural number. The functions 



(1) 



(2) 



are defined making use of the Hermite polynomial iJ„. Writing -ff„(x) in terms of its generating function, _ff„ (x 
^TT exp \2xt — t^] , it is possible to write this state (to use a quantum mechanical terminology) as [8] 

3 <^ i 



/•„(j-) = 1 |1^3 ( A _ li, A ) exp [t^l 



(3) 



t=o 



where 



d^ (z, r) = 2, 6xp [iTrra^] exp [Iiriaz] , Im(T) > 0, 



(4) 



is the Jacobi 793-function, following Vilcnkin's notation [13]. In this notation the basic properties of this even function 
read as 



^3 {z + m + fiT, r) = exp [—iTTTn ] cxp [— 27rmz] ^3 {z, r) 

-,2- 
z?3(2,zt) 



T ""^ " cxp 



TTZ 
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(5) 
(6) 



emphasizing its period 1 and quasi-period r. A beautiful consequence of (6) is that this function can be written as a 
sum of Gaussians 



^Ahi^)^'' ^ '^p 



-7r(y) [aL + zf 



(7) 



a form in which the width — becomes apparent. Property (6) also provides an easy way to obtain the additional 
identity (also given by [6]) 



/n(j) = eHr E °^p 



TT o 27ri 



which is in fact a generalization of Eq.(7) (if one compares it to Eq.(2). 



jot 



Hn (ea) , 



(8) 



III. ORTHOGONALITY OF THE /„'s. 



According to Eq.(l), the functions {fn{j)} are eigenstates of the DFT with associated eigenvalue i". Mehta has 
conjectured that {,fniJ)}n=o i^ ^^ orthogonal set, and thus complete, over a finite set of N points. This reasonable 
conjecture, quite surprisingly indeed, does not hold for arbitrary N (it holds for large N). As, in the following, 
evidence will be collected against the original conjecture, details shall be kept to a level higher than usual. 

Let (/„, fm) denote the inner product 

{fnjrn) = ^.fn*{j)fm{j) 



J=0 



N-1 00 

.2(-z)"+'"E y: cxp 
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H„ (ea) Hm (e/3) 



The sum over {j} is a realization of the modulo N Kroenecker delta. 
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"a,l3 — 



1 a^ 13 (modiV) 
a^ f3 (modiV) 



thus 



(/„, U - 2^(-z)"+" J2 €i exp [-^(a' + /?')] ffn (ea) iJ„ (e/?) , 

a,/3— — 00 



The well known identity. 



exp 



r 



Hkix) 



dy cxp 
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together with the sum over {/?} leads to 

(/n,/m)= ^ / dydz cxp 



a,7— — oo 



'-{y^ + z^) + iyea + ize {a + ^N) 



Hn{y)Hm{z), 



where the infinite sum on {7} covers the equivalence class present in S^^ L Now, the sum over {a} by its turn is a 
realization of a modulo 27r Dirac delta, thus, with the integration over {z} and convenient changes of variables, 



00 />oo 

{fn,fm)^'2TT ^ / dy BXp 



-y 



+ i(yeN — vkN)^ 



Hn[y~—jH,n[y + —j 



7,1; — — 00 

where again an infinite sum is introduced due to the modulo 2tt delta. 

The above expression is rather elucidative. It is not hard to realize that the infinite sums over {7, v} are a direct 
consequence of the equivalence classes brought in by the modulo N Kroeneckcr delta present in Eq.(9). For large 
N, the term corresponding to 7 = u = becomes increasingly important, and a simple check shows that this term 
is exactly 5„„i. Thus, as expected, the limit N ^ 00 recovers the usual harmonic oscillator results. For finite (and 
small) A^, however, all terms in the above summation must be taken into account. 

Following then, the sum on 7 is seen to be a realization of the modulo 27r Dirac delta, Jl^^' (y^N — vnN), and after 
a change of variables one has 



(/n, fm) = 27re ^ exp 



2tt 



Nv 



Nttv^ 
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/^,ij — — 00 

Again, summation over on {/i} must be included to account for the 27r periodicity of the Dirac delta. Splitting the 
sum on v in two sums, over the odd and even integers and shifting the sum on /i by vN results in 

2tt 



{fn,.fm) = 27re ^ exp 



fj..v— — 00 

00 






+2TTe 2_, 6xp 

/^,ij — — 00 



Hn{en - eNv)Hm{en + eNv) 

Hn{en - eNv)Hra{e{n + vN + N)). 



27r, N ^^ , , ,, 

— (Ai+yf-2^iV(« + 1/2)2 



Denoting the second term above by (/„, fm)odd, ii N = 2h + k, where the binary variable k controls the parity of N, 
then 



(/n,/m)odd = 27re ^ exp 



2tt 



{n + h + k/2f - 2itN{v + l/2f 



fl.V—~QO 

and yet again shifting the sum on /i by /i + k/2 and the one on i^ by 1/2, 

"" " 27r 



Hniefi - eNv)Hrnie{fi + vN + N)) 



(/n,/m)odd = 27re ^ (fc) ^ (l)cxp 
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Hn{e{fi-Nv))H^{e{^i + vN)), 



where now X]^-oo(^) denotes a sum over the integers (half-integers) ii k = {k — 1), so that back to the general 



expression, 



27r 9 „ ,^ 9 
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ifnj-m) = 27re ^ cxp 
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OO ( 



i/„(e(Ai-iVt;))i/„(e(/z + i;7V)) 

i/„(e(Ai-iVt;))i/„(e(M + t^iV)). 



27r 9 „ ,T 9 
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Now, recourse to the Hermite polynomial's generating function gives 
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27re 
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E 



cxp 



Qtn Qgra 

V /J,!' — — OO 
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/x^ + 2/ie(t + s) - 2ei/iV(t - s) - f 



2 s^ 



27r 

27r 
1? + 2fie{t + s) - 2ez/Ar(t - s) - t'^ - s"' 



t=s=0 



The sum on {/i} results in a ^3-function in the first term, and a ??3 for /c = or a ??2 for k = I in the second. The 
sum on {ly}, by its turn, gives i?3-function in the first term, and a ??2 in the second, as 



Qn Qn 

ifnjm) = 27re-— -■ 



" dt" 9s™ 
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TT N / \ -K 



exp [-i^ - s^] 
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Using the basic properties, 



^3(z,ir) = T ^ " exp 



92 {z,iT) = T 



- ^-1/2 



exp 
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one gets 



i/rij Jn 



2^3/2 5« ^m 

TV at" as™ 



^^(t + s) 2^\ /e(t-^^\ 
"^ ' TT 'iV/ M 27r '2iV/ 



^ * V I" iVy V 27r ' 2nJ 

Finally, compact expressions can be achieved with 
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and for fc — 1 
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Again, the limit N ^ oo easily recovers the usual results, as the i factor inside the i^-functions guarantees that, in 
this limit, only a term proportional to ■§[k-§-^ 6xp [— 4is]|j^j,^Q survives. Anyhow, with the above expressions any 
term (/„, /„) can be calculated as a sum of ^-function derivatives evaluated at zero. The particular situation m = Q, 
for example, for N even, is quite instructive. In this case 
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and its is immediate to see that all n =odd terms are zero. For n = 2 (and for all even numbers not multipliers of 
4), the symmetry of the binomial term and the multiplicity of the powers of i lead to a pairwise cancelation of all 
non-zero terms. For n = A (and its multipliers), the situation its different. The simplest case is n = 4, 



(/4,/o) = 



(/4,/o) 



r3/2 JU [A 
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tt' N 
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3=0 ^■' 
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This term (with proper normalization) goes to zero quite fast with increasing A^. In fact, for A = 10 it is already of 
order of 10^®. On the other hand, it is immaterial to discuss the case A = 4 (or smaller), as in this situation the 
distinct eigenvalues of the Fourier operator are enough to guarantee orthogonality of the whole set. Considering all 
this, it comes down to, literally, half a dozen different values of the dimensionality A^ (the range [5, 10]) for which a 
significant deviation from the "expected" results (that is, orthogonality) can be observed. 

IV. DFT AND WIDTH INVERSION 

Starting from the own definition of the i^a-function, Eq.(4), with ^ G R, a fractional shift of the "d^ function can be 
calculated. 



k i^2 



^^^'+N^^]= E °^p 



A^^ 



cxp 



' A 



*(^+^.f)=i:i:»p 



exp 



where k is an integer. The sum over {a} can be broken into modulo A^ equivalence classes as 

N 

N' N , 

' j=0/3=-oo 

Conveniently regrouping the terms one gets 

N I ox, 
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A^ 
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exp [2TTi/3 {i^^j + Nz 



X exp 



2m 



'a^ i +'^^iiz + ^ik 



where the term inside the brackets can be identified as ^3-function, 

c2x N 



J=0 



TT 2 2 r. 27ri 



Use of property (6) leads to 



N 



exp 






and taking advantage of the Fourier coefficients exp [^J^] it is easy to obtain the inverse relation 
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(10) 



(11) 



Particular cases of these equations are most interesting, and a lot of peculiar relations can be obtained with the 
different possible choices oi z,k and ^. Two straightforward examples are: First, putting z = in (10), 



N 



k_ i^\ _ _}_^ . (J_ _±_ 



exp 



2m ' 



(12) 



and, according to Eq.(7), the ^3-function on the left hand side has width ^, while the one under the action of the 
DFT has width ^^^. This property is the obvious discrete counterpart of the well known behavior of the Gaussian 
function under the usual Fourier transform. 

The case fc = 0, by its turn, after some manipulation gives 
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^3(iVz,^A^a = ^^;^^3(.+ ^/^' 



TV' N 



A. Application 



With the above results it is possible to generalize the result of [6] in a straightforward way. Introducing 



/-aO-Jf |r^3(^-^af )exp[.- 



t=o 



its DFT can be directly calculated. 
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and use of Eq.(6) together with change of variables from t to it leads to 
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thus 
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/„(A:,r') = *"Ee^P 
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with reproduces Eq.(l) for ^ = 1. From this relation, most identities obtained in [6] may also be generalized. 
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V. TWO VARIABLE'S DFT 

Yet another generalization of the main resuh of [6] regards a two variable DFT, which, for the sake of briefness, 
here it will be merely verified. Apart from the obvious product solution ,fm{j).fn{l), if one considers the quantity 

N-l 



Fm,n{j, = X! ■frn{k)Jn{k - I) CXp 
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which obeys 



[Fm.nU^ 0)* = Fm.nU^ exp 
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use of Eq.(l), and some simple manipulations lead to the non trivial result 
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As in the one variable case, this states obey 

Af-l 

^ |^m,«(j, 01 \Fm',n'{j^)\ =^m,m'&n.ri' m + n ^ to' + n' (mod4), 

which imply a multitude of relations involving derivatives of the ?93-functions (or the Hermite polynomials) . Motivated 
by the previous section, it should be investigated wether this relations holds for m + n = to' + n'(mod4). 

VI. FURTHER RELATIONS INVOLVING THE WIDTH 



So far it has been seen that to break up the infinite sum present in the definition of the Jacobi ^3- function leads to 
interesting properties of this very function. In order to further explore this technique, from Eq.(7) it is straightforward 
to write 

2' 



l?3 



Z'T 



i 



E 



cxp 



a— — 00 



-Ltt 



iF i 



(14) 



with L a positive real number. Choosing ^ integer, it is possible to break the sum over {a} into modulo ^ equivalence 
classes 
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and the infinite sum can be identified as a "d?, 



z i^^ \ 1 ^— i /^ -^ + JF i 
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It is quite interesting to put z = z^ above and observe that 
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which, for the particular case C = 2 gives the well known result 

z i 



^3 


/2z 

u 


Ai 




(2z 


Ai 


1)3 


\ L 


T 



^3 

^3 



L' L 
z i 
L'L 



Similar reasoning would lead to the complementary relation 
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And again, the particular case ^ — 2 gives 
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Equations (15) and (16) can be combined to provide an alternative width inversion relation 
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VII. CONCLUSIONS 

The results here presented seem to argue in favor of one basic point: The Jacobi i?3-function, together with the 
DFT, plays, in finite dimensional spaces, the same role played by the Gaussian function in conjunction with the usual 
Fourier transform. Concerning quantum mechanics, Schwingcr has already noted that, if the families of eigenstates of 
two different observables are connected via DFT, then those observables share a maximum degree of incompatibility 
[5]. In this connection, the width inversion relation obeyed by the /n(j, C) functions strongly suggests that one may 
able to construct, for finite dimensional spaces, states which behavior resemble that of the continuous minimum 
uncertainty states. 

However, such a reasoning meets an important hindrance if one considers that the orthogonality of the DFT 's 
eigenstates ultimately fails. It is a fact, however, that with increasing N it becomes, in a numerical sense, true, and 
in this case the N ^ 00 limit is reached, as witty as it may sound, somewhere near a dozen. This fact may illustrate 
a true finite dimensional idiosyncrasy, or it might lead one to look for the possibility of finding different sets of DFT's 
eigenstates, an issue which is matter of current research. 
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